Magnetic relaxation of type II superconductors 
in a mixed state of entrapped and shielded flux 
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The magnetic relaxation has been investigated in type II superconductors when the initial mag- 
netic state is realized with entrapped and shielded flux (ESF) contemporarily. This flux state is 
produced by an inversion in the magnetic fleld ramp rate due to for example a magnetic field over- 
shoot. The investigation has been faced both numerically and by measuring the magnetic relaxation 
in BSCCO tapes. Numerical computations have been performed in the case of an infinite thick strip 
and of an infinite slab, showing a quickly relaxing magnetization in the first seconds. As verified 
experimentally, the effects of the overshoot cannot be neglected simply by cutting the first 10-100 
seconds in the magnetic relaxation. On the other hand, at very long times, the magnetic states 
relax toward those corresponding to field profiles with only shielded flux or only entrapped flux, 
depending on the amplitude of the field change with respect to the full penetration field of the con- 
sidered superconducting samples. In addition, we have performed numerical simulations in order to 
reproduce the relaxation curves measured on the BSCCO(2223) tapes; this allowed us to interpret 
correctly also the first seconds of the M{t) curves. 

PACS numbers: 74.25.Ha, 74.25.Qt, 74.72.-h 
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I. INTRODUCTION 

In type II superconductors, at temperatures T ^ 0, 
the magnetization relaxes approximately logarithmically 
on time (t) because of the thermally activated motion 
of vortices (flux creep). This behaviour can be under- 
stood, at first sight, within the Anderson Kim model 
(AKM)iiSi^. In conventional superconductor the ex- 
perimental results are well reproduced in the framework 
of the AKM, whereas in high temperature superconduc- 
tors (HTS), deviations from the logarithmic decay are 
observed, especially in Bi-based materials^iSiL^. Several 
models have been proposed in order to explain the non- 
logarithmic relaxation^iiSiiiiiSiH. The theory of collec- 
tive creep, extensively reviewed by Blatter et ali^, pre- 
dicts that the current density (J) relaxes according to the 
so called "interpolation formula" . As in the case of the 
Bean fully penetrated critical state, the magnetization 
can be assumed proportional to the persistent current, 
leading to: 



the equation (Q) immediately leads to: 
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where Mi is the initial value of the magnetization, ks is 
the Boltzmann constant and Uo is the pinning activation 
energy. The exponent /i is a parameter and its value de- 
pends on the different creep regimes; to is a characteristic 
time depending on temperature, magnetic field, sample 
geometry and the fluxon attempt frequency for jumping 
out the pinning centres. By defining the normalized creep 
rate (S): 
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The equation Q is employed to evaluate experimentally 
the pinning activation energy and the exponent fi. For 
this reason, magnetic relaxation measurements are ex- 
tensively used to investigate the flux creep in supercon- 
ductors (for a review, see^ and references therein). 

Usually, in magnetic relaxation measurements (M(t)), 
an external magnetic field Ha ramps up to a fixed value 
Hq with finite sweep rate dHa/dt, then the magnetiza- 
tion is measured as function of time (typically for about 
lO'^ seconds) keeping the external field at the fixed value. 
Ramping the external field up to Hq, that is chosen 
higher than the full penetration field Hp of the super- 
conductor, screening persistent currents (clockwise with 
respect to the external field versus) flow everywhere in 
the superconductor. If the magnetic field is firstly in- 
creased and then slightly reduced, both clockwise and 
counterclockwise persistent currents flow in the sample. 
In this case, the measured magnetization results from 
a region with entrapped flux close to the surface and a 
region with shielded flux in the inner part of the super- 
conductor (ESF state). 

This complicated state can be easily generated when 
the external field ramp is stopped and a magnetic field 
overshoot occurs. This means that, at the nominal stop 
of the external field ramp, the field exceeds the target 
value Hq, reaching it usually after few seconds. This 
overshoot can produce an entrapped flux zone close to 
the surface, which can appreciably affect the relaxation 
process. In particular, Jirsa et al»i^^ showed that, for 
a superconducting slab of thickness 10~^ m in a parallel 
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field Hq = 0.5 T, an overshoot of only 1.5 mT leads to 
an initial magnetization M°^, whose value is about one 
third of the one computed in the absence of the over- 
shoot. However, the depressed magnetization Mov{t) 
relaxes with time converging to the ideal Mid{t) curve 
computed in absence of overshoot. Therefore, the initial 
value of the magnetization, occurring in the absence of 
the overshoot, is determined approximately by extrapo- 
lating it from the long time M{t) curve. 

However, starting from the ESF state, the field profile 
evolution that leads to the joint of the two curves is still 
unclear. On the other hand, it is not possible to deter- 
mine experimentally when the M{t) curve approaches to 
the ideal relaxation and, thus, it is usually adopted the 
experimental procedure of cutting the first 10-100 sec- 
onds in the experimental M{t). 

In order to justify this experimental procedure, we can 
consider a slab of thickness 2 d and critical current den- 
sity Jc analyzed in the framework of the Bean model. If 
an overshoot occurs after the application of an external 
field higher than the full penetration field {Hp = Jed), 
the magnetization of the slab in the framework of the 
Bean model, is: 

M = M,„ + M,n (4) 
Me„ = il/i)iHl/Hp) (5) 
M,, = -{l/2){{Hl-Hl)/Hp) (6) 

where M^n is the magnetization due to the entrapped 
fiux, Msh is the magnetization due to the shielded flux 
and Hov is the amplitude of field overshoot. If Hov ^ Hp, 
the magnetization due to the entrapped flux is small and 
thus it can be considered negligible after a long enough 
time. In a low Tc superconducting slab, with d = 0.1 mm 
and Jc = lO^^A/m^, the full penetration fleld is Hp = 
0.63 T and the usual characteristic time Iq is about 10 
seconds. Therefore, for a few mT overshoot, it is com- 
monly believed that the experimental M{t) measured 100 
seconds after the nominal stop of the external magnetic 
field resembles the relaxation from a fully shielded state 
(or a fully entrapped state). Nevertheless, depending on 
the temperature and the applied magnetic field. Hp can 
become comparable with Hov, drastically affecting also 
the long-time magnetic relaxation. 

To extend the relaxation analysis to the time window 
affected by the overshoot, Jirsa et alii^ii^ have shown 
that it is possible to use magnetic hysteresis loop data 
measured at different field sweep rates. They have shown 
how the magnetization measured at different sweep rates 
can be converted into magnetic relaxation data, substan- 
tially extending the time window to the short times, typ- 
ically down to 10~^ s. 

Other complications in the analysis of relaxation mea- 
surements can also arise from the sample geometry and 
the anisotropic properties of the material. In fact, in HTS 
samples, magnetic relaxations are usually measured with 
the field orientation perpendicular to the largest face of 
the sample. In this geometry, the demagnetization effects 



could be neglected only for measurements performed at 
fields much higher than Hp. Since an overshoot changes 
the direction of the current and the magnetic field value 
on the edge of a flat superconductor, geometry effects are 
supposed to be significatively altered in the magnetic re- 
laxation measurement. 

In this work we have investigated the magnetic relax- 
ation starting from a state with entrapped and shielded 
fiux. 

In the next section, we will discuss the integro- 
differential equation employed in the numerical compu- 
tation of the M{t) curves. In the Sect ion ITTll we show the 
numerical simulations of the magnetic relaxation and the 
time evolution of the field profiles for samples in shape 
of slab and thick strip. 

The magnetic relaxations in BSCCO(2223) have been 
experimentally investigated when the effects of a mag- 
netic field overshoot in the M{t) are not negligible. Fi- 
nally, in the Section llVI the experimental measurements 
are analyzed and compared with the numerically com- 
puted results. 



II. NUMERICAL COMPUTATIONS 

In order to analyze the magnetic relaxation of a su- 
perconductor in an external magnetic field Hq, we nu- 
merically solved an integro-differential equation for the 
current density J in a slab in parallel field and in a thick 
strip in perpendicular fieldAS. As developed by Brandt 
in a series of worksiSiiLiSii2iS£, in a long strip of width 
2a (along y axis) and thickness 2d (along z axis) placed 
into a homogeneous magnetic field, perpendicular to the 
largest face of the strip, the applied field induces surface 
and bulk currents. The current flows along the sample 
length (i.e. x axis) due to the symmetry of the strip. 
The induced current density J — J{y, z)'i generates a 
magnetic field H which has y and z components. In this 
model it is assumed that B = /^oH and thus, Hd and 
the reversible magnetization {M^ev) are neglected. Since 
B = V X A, where A is the vector potential, it is possible 
to write for this geometry a 2D Poisson equation in the 
Coulomb gauge 

/^oJ^-V^A (7) 

The current density flows only in the strip and thus the 
vector potential could be written as a sum of two terms 
A = Aa + Aj, where Aa is the vector potential related to 
the applied magnetic field, [Aa = [r x B]^ = yBa), and 
Aj is related to the current induced in the strip. Since 
Ba is constant in the specimen, the general solution of 
the lO is: 

^(r) ^ -^^o I d\'Q[v, v')J{v' , t) ~ yBa (8) 
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where r = {x,y), r' = {x',y') and Q{r,r') is the integral 
kernel defined as: 



Q(r,r') = :^log 



r — r 



(9) 



in which rg is an arbitrary constant length that can be 
chosen equal to 1. The integration is performed on the 
cross section of the strip S. The current density, is ob- 
tained formally fromiSi: 

J(r,t) = -— / d\'Q~'{r,r') U(r',t) -ysj . (10) 
Mo Js' L J 

Here Q^^{r, r') is the inverse kernel defined by: 

/ dVQ-i(r,r')g(r',r") = <5(r-r"). (11) 
Js 

By using the relation E = —Wx<t> — A where (j> is the 
scalar potential, we obtain 

j{r,t) = -f d\'Q-\r,r')[E{J)-y'Ba{t)] (12) 
Mo Js' 

In the limit a (slab geometry) , the previous equation 
becomes an one-dimensional equation: 

j(r, t)^ — I dy'Q2,{y, y') [E{J) - y'Ba{t)] (13) 
Mo Jo 

Taking into account the symmetry of the strip and slab 
geometries, the kernel in the case of the strips is given 
by: 



_ 1 iyl+zl){yl+zX) 



where y± ^ y±y' and z± — z±z'. For the slab it results 

Qsiab = ^(ly ~ y'\ - \y + y'\) = - "^Hv, v') (i5) 

In our simulations, we do not consider a transport current 
but only an external magnetic field and for this reason 
the term y x4' has been dropped out. To solve the integral 
equation for J we choose the widely used relationSi: 



Jc 



(16) 



where is the critical current density. However, 
the Brandt method can be used with different E — J 
relationship^^. 

The current density profiles in the strip has been ob- 
tained by integrating the equation (|12|l . whereas for 
the slab the equation (|13|l has been solved. For the 
strip, the functions J, E have been tabulated on a 2D 
grid with equidistant points yk ^ {k — l/2)a/Ny (k = 
-Ny + 1, 0, ■■■,Ny) and zi = {I - l/2)d/N^ {I ^ 
-N^ -I- 1, • • ■ , 0, • • ■ , iV^), where = d/aNy is cho- 
sen. Labelling the points {yk,zi) by an index i, with 



I = 1,2, • • • ,iV and iV = NyN^, the function J{y,z,t) 
becomes the time dependent vector Ji{t) with N co- 
ordinates and E{y,z,t) = Ec{J/Jc)^ becomes a vec- 
tor with N coordinates. Moreover, the integral kernel 
Q{yj z^ y', z') becomes an x matrix Qij . 
The numerical form of the equation Hll)|l is 



J»(i + At) 



At 



N 



for i = 1 , • ■ 



E,it)-yB 

(17) 



where Ay = a/Ny and Az = d/N^ are respectively the 
steps in the 2D grid used to tabulate the cross section 
of the thick strip. The numerical integration of the ID 
equation for a slab follows similar rules. 

The time integration of this system of non-linear dif- 
ferential equations for Ji{t) has to follow some prescrip- 
tions. First of all, the integration starts with the initial 
condition Ji(0) = 0; in addition, the time step At is cho- 
sen inversely proportional to the maximum value of the 
resistivity pi = Ei/Ji. Brandt uses the following rela- 
tion in his computations: At = ci/[max{pi{t)) + C2] with 
ci — 0.3/ {Nyu), n is the exponent in the E — J law and 
C2 = 0.01. In our computations we do not use normalized 
quantity and we have observed that this choice depends 
on the value of Jc and the time derivative of the external 
magnetic field. In our computations we used different 
values for ci and C2 in order to make stable the numeri- 

Bn where Br, is the 



cal algorithm: a = 0.003/ [(iV» 
temporal mean value of Ba and C2 = 1. 

Finally, Qij — In jr^ — r^ j has a logarithm divergence 
when Ti approaches Vj . In order to avoid this singularity 
for i = j the expression for the kernel is changed with 



(l/2)ln[(r, 



e^] wherai^: 



= exp[ln(A?;2 + Az^) - ln(4) - 3 + 

Ay ,Az. .Az. .Ay,, 

+ arctani-— ) -f (- — larctani— — . 
Az ^Ay' ^Ay' ^ Az'^ 



In our computations, the magnetization is calculated by 

M = A dy dzJ(y, z)y for a strip (18) 
Jo Jq 

M = 2 / dyj{y)y for a slab (19) 



Since the magnetic relaxation is simulated on 10'"' ^ 10^ 
seconds, we reduce the number of computed points cal- 
culating the {ti, Mi) data accordingly to the relation: 

t, = t,_i +exp(log(tfl)/7Vfl) (20) 
M, = M{t,) (21) 

where tj^ is the total time of the computed relaxation and 
Nji is the total number of the computed data. 
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FIG. 1: Time dependence of the external magnetic field 
during the field overshoot. The time origin corresponds 
to the nominal field stop. A triangular overshoot (dashed- 
dotted line) and an overshoot given by the function AHov ~ 
-H'ou(t/iot)m)^ exp(2(f — t/tovm)) (dotted line) are shown. In 
the inset, the time derivatives of the two overshoot functions 
are shown. 



III. NUMERICAL RESULTS 

In this section we discuss the numerical results ob- 
tained for the slab and the thick strip. In our computa- 
tions we have used a strip with aspect ratio (a/d) equal 
to 10 and 2a = 10~^ m and a slab with 2a = 10"'* m, 
with the critical current density (J^) ranging from 10^ 
A/m^ to 10^ Kjvr? . The current-voltage characteristic 
is the usual power law given by = Ec{J/ JcY^ where 
Ec ~ V/m and the employed exponent n is chosen 
equal to 15 for the large creep case and n = 105 in the 
Bean limit case. 

In order to study the relaxation from a ESF state, dif- 
ferent magnetic field ramps have been taken in account. 
For each ramp, the external magnetic field Ha increases 
linearly on time, with a sweep rate {Ha) of 1 mT/s, up 
to a nominal fixed value Hq. The time when Ha has 
nominally reached Hq is taken as time origin of the mag- 
netic relaxation. As the external magnetic field reaches 
Hq different situations are taken into account: 

a) Ha is stopped immediately (ideal case); 

b) Ha has a triangle overshoot (triangle overshoot); 

c) Ha has an overshoot with a smoothed field stop 
(exponential overshoot) . 

In the case b), the magnetic field increases in tovm sec- 
onds by an amplitude Hov, then it decreases by the 
same quantity in the subsequent tov seconds (triangle 
overshoot). After this, the external field is immediately 
stopped and the magnetic relaxation starts. In the case 
c), the overshoot has been simulated by means of the 
function Fov{t) = Hovit/tovmY Gxp{c{l - t/tovm)); for 



FIG. 2: Magnetic relaxation curves computed for different 
magnetic field ramps in a thick strip. 



t — tovm the overshoot reaches the maximum value. The 
two different functions employed to simulate an overshoot 
are shown in Fig. ^a,). For the triangular overshoot, we 
set Hov — 1 mT, tovm = 1 s, tov = 5s. In the case of 
exponential overshoot, we used Hov — 1 mT, tgvm = 1 s, 
c = 2. In the inset of the same figure, the time deriva- 
tive of the overshoot functions are plotted, since the field 
ramp derivative is actually used in the integration of the 
diffusion equation. 

We have initially computed the magnetic relaxations 
for a strip in perpendicular magnetic field (perpendic- 
ular geometry) by simulating a case analogous to the 
one discussed in the work of Jirsa et aU^. In our com- 
putation Jc = 10^ A/m^ and the critical exponent is 
n = 15. We are considering a superconductor with large 
critical current density but with large creep. The ex- 
ternal field ramps with a sweep rate of 1 mT/s up to 
0.2 T, which is a value well above the full penetration 
field of the strip. Indeed, looking at the field profile we 
have verified that the strip is fully penetrated for fields 
higher than 0.10 T. As shown in Fig. |21 also if the over- 
shoot does not occur in the field ramp, the magnetiza- 
tion decays non-logarithmically, especially at short time 
(< 10 s). This result is expected due to the power law 
in the E — J relationship which involves a logarithmic 
dependence of the pinning energy on the current density. 
In the same figure, a magnetic relaxation curve is shown 
as computed for a field ramp which has a triangular over- 
shoot. In this case, the external magnetic field ramps up 
to 0.2 T. After this, the field overshoot occurs with an 
amplitude of Hov = 1 niT. The field overshoot reaches 
its maximum one second after the external field should 
have been stopped at the nominal target value. The ex- 
ternal field goes down to the nominal value of 0.2 T after 
5 seconds. Also in this case, when the magnetic field is 
stopped at the fixed value of 0.2 T the time derivative 
of Ha is instantaneously zero. A more realistic situation 
have been considered by computing the magnetic relax- 
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ation for the case c) where Hov = 1 mT, tovm = 1 s and 
c = 2. The case c) is effectively reahzed in experiments, 
where the field cannot be stopped instantaneously and 
the overshoot shape is rounded. 

As shown in Fig. [2 the magnetization curves in the 
cases b) and c)have an initial values Mi larger than in 
the ideal case. In fact, when the overshoot occurs, the 
magnetization does not relax during the first seconds, 
since the magnetic field continues to increase. The largest 
value of the initial magnetization is obtained in the case 
of an exponential overshoot; indeed, the electrical field in- 
duced in the superconductor in the first seconds is larger 
than in the other cases (see also Fig. 0. When the 
external magnetic field rate reverses, the magnetization 
quickly decreases, because of the flux coming out from 
the surface, and after 5 seconds M has lost the 12% of 
the initial value. The decay during the first 5 seconds 
depends on the shape of the field overshoot as function 
of the time. In the triangular case, the magnetization 
curve shows a convex concavity, whereas in the case c) 
the curvature is concave. After 5 seconds, the external 
field is practically constant and the magnetic relaxation 
effectively starts; for t larger than 100 seconds the three 
curves join together. These computations confirm also 
in perpendicular^eometry, the results found for parallel 
geometry in Ref;;14. However in this case the field over- 
shoot amplitude is 1% of the full penetration field. In the 
next section we will consider situations where the induced 
ESF state strongly affects the magnetic relaxation. 
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FIG. 4; Time evolution of the magnetic field profiles in a slab; 
the initial field profile is achieved without field overshoot. 



1.00 



QQ 
CD 



0.99 









fs' ■ 










2 s 








"'A 
TO? 


3 s 




























Is- 










2 s : 






VP' 




3 s ; 






'tl : 1.0C 


Os 










0.5C 


y/a 



-1.0 



-0.5 



0.0 0.5 

y/a 



1.0 



1.5 



T — I — I — I — I 




Llllll I I I I I I i 

10' 10° 10' 10' 10' 10' 10' 10' 



t(s) 

FIG. 3: Magnetic relaxation curves computed for different 
magnetic field ramps in a slab. The dotted line is the relax- 
ation given by the analytical formula reported in the text. 

Here, we discuss the magnetic relaxation starting from 
an ESF state in the case of a slab in parallel field. In 
Fig. 131 two computed M(i) curves are shown; the initial 
magnetic state is obtained by ramping the external field 



FIG. 5: Time evolution of the magnetic field profiles in a 
slab, when a field overshoot occurs. The dotted lines repre- 
sent the fiux profiles which fully resemble the ideal profiles. 
The dashed lines represent the profile during field ramp rate 
reversing. Continuous lines show the field profiles in the time 
windows where the magnetization is nearly constant. In the 
inset a detail of the profile close to the slab surface is shown. 

both in the ideal way (without overshoot) and with an 
overshoot of 1 mT (dashed curve) In the same figure, 
it is shown the magnetic relaxation (dotted line) for a 
superconducting slab, according to the relation given in 
Ref. [a where it is assumed that the pinning energy 
depends logarithmically on the current density: 

M(t) = M(0)exp ^-iln ^^^^ (22) 

The dimension of the slab used for the computations is 
2a = IO~'*m, the critical current density Jc = 10^ A/m^ 
and the exponent n = 15. In this case the full penetra- 
tion field of the slab is ~ 6.3 mT and thus it is of the 
same order of magnitude with respect to the overshoot 
(1 mT). As shown in Fig. |31 the computed ideal curve is 
approximated quite well by the analytical relation in the 
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time range from 10 to 10^ seconds, whereas it wanders 
off each other at very short and very long times. On the 
other hand, we observe as the overshoot has effects on 
long time up to 10^ s (dashed cm've). In the first 5 sec- 
onds, the magnetization looses 60% of the initial value 
due to the inversion of the flux profile close to the slab 
surface. In the subsequent 10^ seconds the magnetiza- 
tion practically does not relax, and after this time the 
relaxation rate increases. After 10^ seconds the mag- 
netization computed with an ideal ramp and the curve 
computed with a field overshoot take the same value. 

At this point, it is necessary to investigate if the mag- 
netization computed for time larger than 3.0 10^ s in 
both the cases, corresponds to the same magnetic state. 
In order to answer this question, we have computed the 
magnetic field profiles as a function of time. In Fig. ^ 
and in Fig. [S] the field profiles computed for both the 
cases are shown. In particular, in Fig. ^ the profiles of 
the relaxation in a slab are shown, reproducing the usual 
Bean results. On the other hand, the profiles computed 
in the case of a relaxation from an ESF state, obtained by 
using the exponential overshoot, (Fig. [Slshow that dur- 
ing the first 5 seconds the profile changes (dashed line) 
as a consequence of the field decreasing. The evolution of 
the profiles during the first 5 seconds has some difference 
in comparison with the classic Bean profile, where Jc is 
constant and independent on the applied electrical field. 
In our case, while the flux is expelled on the surface, in 
the inner part of the slab the profile relaxes. This occurs 
because of the finite exponent n which leads to a large 
creep. On the contrary, for the Bean model, the field pro- 
file, in the inner part of the slab, remains frozen during 
the field decreasing. 

Starting from the fifth second the field profile relaxes 
overall in the slab and after 10^ seconds the magnetic 
profile becomes the ideal one. In Fig. |5lthe field profiles 
which resemble the ideal ones are shown by dotted line. 
By means of our numerical simulations we have shown 
that the same magnetization value found in the two M{t) 
curves corresponds to the same magnetic state. In Fig. 
|31 we observe also that the maximum of the field profile, 
due to the field ramp rate reversing, moves towards the 
slab edges during the relaxation. At the same time, the 
entrapped magnetization is reduced down to zero. There- 
fore the ESF state has relaxed towards a fully shielded 
state. 

Increasing the amplitude of the overshoot, we expect 
that the ideal relaxation and the relaxation from a ESF 
state will coincide at longer times. Nevertheless, as the 
region with entrapped flux prevails on the shielded re- 
gion, the flux profile relaxes towards a fully entrapped 
state. 



B. ESF state in strip 

In order to analyze the effect of the sample geometry 
on the relaxation, we considered the case of a strip in 
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FIG. 6; Evolution field on time for a thick strip. The field 
profiles computed in the ideal case are shown on the upper 
frame. The field profiles computed when an overshoot occurs 
in the field ramp are shown on the lower frame. 





STRIP / 








STRIP 


1.01 


^ = // 




1.010 


5 s/ ^^ 








\\\ Os 


















lO^^S'' \ 




m° 
























Id's " .-'v 












CO 1 005 
















io's,-'',,--'' 
















1.00 








lV s 





0.5 0.5 

y/a y'a 



FIG. 7: zoom of the previous figures, in the region close to 
the strip edge. 



perpendicular field for which the main effect of the over- 
shoot arises on the surface, where the demagnetizing field 
is more intense. In Fig. Elthe magnetic field profiles for 
a thick strip (2a —1 mm, 2b =0.1 mm) are reported; a 
critical current density of 10^ A/ui^ and an E — J ex- 
ponent n = 15 are set. In the upper part of the figure 
we can see the field profile relaxations in the ideal case. 
We can observe that the demagnetizing field relaxes to- 
wards lower magnetic fields on the surface. At the same 
time, the field increases in the inner region and there is 
a boundary, known as the neutral line, where the field 
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evolution of the current density is shown . During the 




relaxation the current changes sign and after long time 
the current distribution in the cross section of the strip 
rebuilds the distribution of a full shielded state. Except 
for the time evolution of the magnetic field on the border 
of the strip, in the perpendicular geometry there are not 
substantial differences respect to the parallel geometry. 
In fact, our computations have shown that in the per- 
pendicular geometry, for Ha > Hp, the demagnetizing 
effects do not affect the time evolution of the magnetic 
relaxation. 



IV. EXPERIMENTAL RESULTS AND 
DISCUSSION 



FIG. 8: Magnetic relaxation curves computed from different 
magnetic field ramps in a thick strip. 
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FIG. 10: Magnetic field ramp with a sweep rate of 3.3 mT/s in 
the time window where a field overshoot occurs, as measured 
by a hall-probe (square and line), and Ha{t) employed in our 
computation (dotted fine). 



FIG. 9: Current density distribution at different times for 
a strip. The current density is computed in the case of an 
overshoot in the field ramp. 



remains constant; it divides the region with entrapped 
flux from the one with shielded flux. If an overshoot of 
a 1 mT occurs the flux, as expected, is strongly reduced 
on the strip edge and the field maximum is located inside 
the strip. In the next 10^ seconds the maximum relaxes 
and moves towards the strip edge where, at the same 
time, the field increases. On the contrary, in the ideal 
case the field on the border always decreases during the 
relaxation. When the maximum reaches the edge, the 
field profile in the strip fully resembles the profile com- 
puted in the ideal case and the relaxation continues as in 
the ideal case. Also in this case, as shown in the mag- 
netization curves in Fig. (SJ the M{t) with and without 
overshoot join together at long times. Also in the perpen- 
dicular geometry the evolution of the magnetic state is 
directed to rebuild a shielded state. In Fig. |51 the time 



Magnetic relaxation measurements have been per- 
formed by means of a Vibrating Sample Magnetometer 
(VSM) equipped with a 16 T magnet. The external mag- 
netic field can be ramped with a maximum sweep rate of 
7 mT/s. When the field is nominally stopped the mag- 
netic field has an overshoot of around 1-^5 mT depending 
on the sweep rate used for ramping the field and this un- 
wanted feature has been used to induce a ESF state in 
our samples. We used a hall probe to measure the time 
dependence of the external field and in the inset of the 
Fig. mithc measured overshoot for our magnet is shown. 

In order to check of validity of our numerical results, we 
have measured the magnetic relaxation on monofilamen- 
tary BSCCO(2223)/Ag tapes prepared by the standard 
PIT technique. We have chosen this kind of samples be- 
cause they allow us to study bulk rectangular samples 
with full penetration fields which can be in the order of 
10 mT even at the lowest temperature i.e. 4.2 K. The di- 
mensions of the superconducting region in the measured 
sample are 3.02 x 0.14 x 4.6 mm'^ and the estimated 
critical current density ranges from 10'^ to 10^ A/m^, de- 
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FIG. 11: Magnetic relaxations measured at different temper- 
atures for a magnetic field -ffo = 1 T (points) and computed 
curves (continuous lines). 



FIG. 12: Field profiles computed for the relaxation at 45 K. 
The profile in the direction of the arrow are computed at 
t =100, 1000, 5000 s. 



TABLE I: critical current densities and exponent n used for 
the fit of the experimental curves. 

T (K) Jc (A/m^) n 



4.2 2.40 10** 20 

15 1.15 10* 19 

25 9.70 lO'^ 13 

35 8.70 lO'^ 9 

45 4.00 lO'^ 8 



pending on the temperature. In this way we can study 
experimentally the overshoot effects as Hp decreases. 

M{t) measurements have been performed with the field 
perpendicular to the sample surface {H \\ c-axis) in the 

4.2 - 45 K temperature range, cooling the sample in zero- 
field (ZFC) for each temperature. The initial magnetic 
state is obtained by increasing Ha with a sweep rate of 

3.3 mT/s, up to 2 T. After this, the field is decreased 
with the same sweep rate down to a measuring field 
IIqHq = 1 T. The field variation of 1 T is chosen to be, 
for any measuring temperature, well above i?p, which is 
evaluated by taking the value of the field corresponding 
to the maximum (in absolute value) in the virgin mag- 
netization curves at 4.2 K. In this way, in absence of a 
field overshoot, a full critical state, with entrapped flux, 
is realized in the superconductorsi^. As the final field Hq 
is nominally achieved, the M{t) data are acquired each 
second for 5000 seconds. 

The M{t), normalized at the initial magnetization value 
^(0), measured at different temperatures, are shown in 
Fig. ^2 In all the curves, a large drop in the magneti- 
zation occurs during the first 11 seconds and this time 
corresponds to time interval during which the external 
field has an overshoot. The behaviour of the magneti- 
zation in the subsequent 5000 seconds depends on the 
value of the temperature. At 4.2 K the magnetization 



decreases slightly, but the relaxation after 5000 seconds 
does not exhibit the behaviour expected for a fully en- 
trapped state. At 15 K and 25 K, the magnetization 
remains nearly constant, whereas at 35 K and 45 K the 
magnetization first takes negative values and then in- 
creases on time. The effect of the overshoot increases as 
the full penetration field decreases with the temperature. 
These measurements show that the magnetic relaxation 
can be still affected by the field overshoot after at list 
100 s. The negative values measured in the M{t) at 45 K 
mean that the shielded flux region in the sample is larger 
than the entrapped one, although the initial condition 
was a fully entrapped state. 

In order to reproduce our experimental results, we 
have computed the magnetic relaxation for a supercon- 
ducting strip with the cross section of our sample. In 
the computations, the field ramp reproduces the experi- 
mental field ramp, with a sweep rate of 0.0033 T/s. The 
overshoot has been simulated by using the exponential 
function discussed in the Section UTTl As shown in Fig. 
1101 this function reproduces quite well the experimen- 
tal overshoot with Hoy = 0.029 niT, tovm — ^-"^ s and 
c = 1.3. In our computation we have to set both n and 
Jc- The exponent n has been evaluated by measuring 
the hysteresis loop at different sweep rate. Taking the 
M values measured at 1 T for different sweep rate (Ba), 
n is given by means of a linear fit of log(-Ba) as function 
of log(Af); the n values reported in Tab. ^ have been 
rounded to the nearest integer. On the other hand, the 
critical current density is a free parameter chosen in or- 
der to obtain the best fit. From our computations, it 
results Jc = 2.4 10^ A/m^ at 4.2 K and 4.0 10^ A/m^ at 
45 K. As shown in Fig. ^2 the numerical computations 
reproduce well the experimental behaviour. In Fig. 1121 
the profile computed at T=45 K are shown. In particu- 
lar, at t=10 s, when Ha is practically constant, it results 
that the magnetic state in the superconductor has both 
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the regions with entrapped and shielded flux. In the next 
5000 seconds, the profile relaxes toward a shielded state, 
which is practically fulfilled at t=5000 s, when the simu- 
lation is stopped. 

Our work shows that the first seconds of the relaxation 
have to be analyzed very carefully in order to estimate 
correctly the creep rate and, thus, extract information 
about the pinning properties of the sample. In fact, our 
results show that it is not appropriate just to cut the first 
seconds of the relaxation curves and extract information 
from the remanent data if the presence of an overshoot 
in the magnet has not been previously considered. 

V. CONCLUSION 

In this work, we have studied the magnetic relaxation 
from a state with shielded and entrapped fiux, generated 
by a field overshoot after the nominal stop of the exter- 
nal field. The magnetic relaxations have been computed 



in parallel and perpendicular geometry. The computed 
magnetization shows a large drop in the first seconds due 
to the flux expulsion from the samples boundary. Af- 
ter long time, the M(t) curves computed with and with- 
out field overshoot (having, thus, as initial condition an 
ESF and a full shielded or entrapped flux state, respec- 
tively) join together. Moreover, our simulations show 
that, during the relaxation, the same value of the mag- 
netization corresponds to the same magnetic state. In 
addition, the experimental relaxation curves, measured 
on BSCCO(2223) tapes, are well reproduced by our nu- 
merical computations, allowing us to correctly analyze 
the M(t) from the instant when the external field is nom- 
inally stopped. 
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